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Abstract 

The Yang-Mills magnetofluid unification is constructed using lagrangian 
approach by imposing certain gauge symmetry to the matter inside the fluid. 
The model provides a general description for relativistic fluid interacting with 
Abelian or non-Abelian gauge field. The differences with the hybrid mag- 
netofluid model are discussed, and few physical consequences of this formalism 
are worked out. 
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Some experimental discoveries and studies suggest that the deconfined quark 
gluon matter is behaving more like a quark gluon plasma (QGP) liquid [1, 2, 3]. 
This fact motivates tremendous works in constructing the non-Abelian fluid models 
like magnetohydrodynamics [4, 5, 6, 7, 8, 9, 10, 11, 12]. 

In some recent models [9, 10, 11], the relativistic hot fluid was described in terms 
of hybrid magnetofluid field which unifies the electromagnetic and fluid fields. The 
unification is represented by the effective field strength tensor, M^,^ = F^^ + m/qS^i, 
combining appropriately weighted electromagnetic and fluid fields. The model has 



* Email 
Email 
^ Email 



lyman@tisda.org, sulaiman@teori.fisika.lipi. go. id 
tpdjun@teori.fisika.lipi. go. id 

handoko@teori.fisika.Iipi. go. id, laksana.tri.handoko@lipi.go.id 



further been generalized to the non-Abehan case [10]. Also its topological structures 
has been investigated in detail to model the glueballs [11]. 

In this Letter, we follow the same line but with different starting point using 
lagrangian approach and imposing certain gauge symmetry on it. This approach 
is intended to improve the Yang-Mills magnetofluid unification model [10] in two 
aspects : 

• Prom field theory point of view, the fluid and another gauge flelds should 
physically be different fields, and then should have independent kinetic terms 
in the lagrangian. Deploying the above effective field tensor would lead to the 
unphysical two point interaction connecting two different fields from the gauge 
invariance kinetic term M^^M'^'' in the lagrangian. 

• The hybrid strength tensor M^i, suggests that the fiuid and another Abelian 
or non- Abelian gauge fields are gauge invariance under a single gauge trans- 
formation and both should belong to the same Lie group. This requirement 
is actually necessary since the formalism relies on the Lorentz force equation. 
This also avoids us to deal with, for instance, the case of non-Abelian fluid 
interacting with an (Abelian) electromagnetic field. 

In order to overcome such problems, we propose local and independent gauge 
symmetries for each field of fiuid and another gauge boson interacting with it, in- 
spired by the unified theory known in particle physics. The model is able to deal 
with some scenarios like : (i) non-Abelian fiuid interacting with non-Abelian gauge 
field under G{n)-p ® G(n)G symmetry; (ii) non-Abelian fluid interacting with an 
(Abehan) electromagnetic field under G{n)-p ® U(1)g symmetry; and (iii) Abehan 
fiuid interacting with an electromagnetic field under U(1)f ® U(1)g symmetry, i.e. 
the original magnetofiuid model. Here F and G denote the fiuid and another gauge 
fields respectively. Note that G(n) could be any n dimensional Lie group like SU(n). 

Following the basic procedure in lagrangian approach, we start from the following 
matter lagrangians, 

r _ / (a„$)"^(c"'$) + ^m|$^$ + for boson 

''-matter — \ 2 ■ 

y i^^^ — m^^^ for fermion 

V{^) is the potential, for example in a typical theory, V{^) = iA($*$)2, and 
^ = 7^9'* with 7^'s are the Dirac matrices. 

Under a local non-Abelian G{n) gauge transformation U = exp[—iTa9°'{x)] ^ 1 — 

iTa9°-{x) with 9^ -C 1, the matter field is transformed as $ = exp[-ira^"(a;)] $ 

or ^ ^' = cxp[— zTa^^"(.T)] ^ with the matter field is in general an n x 1 multi- 
plet containing n elements for n dimension Lie groups as SU(n), 0(n + 1), etc. T"'s 
are generators belong to those Lie groups and satisfy certain commutation relation 
[T", T^] — if'^'^T'^ with /"^"^ is the anti-symmetric structure constant. Anyway, the 
number of generators, and also gauge bosons, is determined by the dimension of 
group under consideration. For an SU(n) group one has ri^ — 1 generators and the 
index a runs over 1, 2, • • • , — 1. The U(l) is realized by a phase transformation. 
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Ta6"'{x) — >■ 6{x). The gauge invariance is then revealed by introducing gauge field 

A^'s which are transformed as A"^ A"^' = + (l/i/)(<9;.^") + /"''^^^A^^, and 
replacing the derivative with the covariant one, = d^ — igT"-A1^, where g is the 
gauge "charge" . Further the gauge invariance kinetic term for gauge boson ^4^ takes 
the form of F^^F""'"', with strength tensor F^^ = df^A^ - d^A^ + gf^^^A^^Al. 

Now, let us consider the most general case of non-Abclian fluid interacting with 
any non-Abelian gauge flclds represented by Giri)^ ® G{n)Q symmetry. Hence we 
can deploy the above gauge principle for each field independently, that is the gauge 
transformations are performed separately with independent phase parameters ^(x)'s. 
Imposing the symmetries, and assigning and ^4^ as the gauge bosons associated 
with fiuid and another gauge fields belong to different group spaces, the covariant 
derivative becomes, 

V,^d, + ig^ + ig^ T^Al , (2) 

where g^ is the "charge" for fluid, while gQ is the gauge charge. 

Finally the total lagrangian with such gauge symmetries becomes, 

^ — -^matter ~l~ -^gauge ~l~ -^int. i (3) 

where, 

1 1 

''-gauge — ' ^ ' 

Ant. = -g^J^^U^^-gGJG,A^' + J:}ir\ (5) 

and S'^j^ = d^V^ — d,yU^ + gpf'^^'^Uj^UI;. Meanwhile, in the case of bosonic matter 
there are additional interference terms in Eq. (5) coming from the gauge invariance 
kinetic term 

£bo«°"^ = gl {^^T^T^^) U^U^^ + gl {^''T^T^^) A^A^^ 

+gFgG {T^T^ + T^T^) U^A'"" , (6) 

The 4-vector "current" is, 

ja ^( -t[id^^yT^^-^^T^{d^^)] for boson 

I *T^7^* forfermion ' ^ ^ 

with X : F, G and = \E'^7o. 

Having the total lagrangian at hand, we can investigate further the dynamics 
of fluid and the contribution of its interactions with another gauge fields. First of 
all, let us derive the equation of motion (EOM) of fluid that is our main interest in 
this Letter. It can be obtained from the lagrangian under consideration through the 
Euler-Lagrange equation in term of U^, 

'"^-^.^ = 0, (8) 
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Substituting Eq. (3) into Eq. (8) yields the EOM, 

(9) 



Sfj^y — QyJy^ for Abelian , 



= qyJ^^ for non - Abelian , 

with the covariant current is, 

=\ -^[{^^,^)^T^^-^^m^^M for boson 
~ I Jp.. forfermion ' ^ ' 



Moreover, the antisymmetric strength tensor in Eq. (9) imphes that d'^J^^ = 0, 
that is „ is a conserved current. 

Following the previous works in [9, 10, 11], we can take fluid fleld to have a 
relativistic form as follow, 

[/; = (US, = ul<t> with u, = 7"(1, -v^ , (11) 

where is relativistic velocity and the relativistic factor 7 = (1 — |vp)~^^^, while 
V is the spatial velocity. On the other hand is a dimension one auxihary field to 
keep correct dimension and should represent the fluid distribution in the system. 
We should remark that the index a in is only to label each flow fleld and not 
summing up all of them. Therefore, in the present model the fluid is described in 
terms of its kinematics and distribution function separately. It captures a compli- 
cated dynamics modeled as a fluid (hke) system. In the case of non-Abelian fluid, 
we consider several flow flelds labeled by the index a. We should remark that, Uf^ 
in Eq. (11) differs slightly with the conventional notation. The reason of taking 
negative sign in the spatial component of will be clarifled soon. However, more 
importantly the form still keeps Lorentz invariance. 

Now we are ready to investigate in detail the EOM for general non-Abehan case 
in Eq. (9). For the sake of simplicity and concerning our current interest, Eq. (9) 
is rewritten in terms of double derivative, covariant current and the rest denoted 
by "total force" = (Fq , F") = F^{U,A). Hence summing up all components we 
simply have. 



d' {d,us - aot/;) - {d,ut - d,u;) = ( + f;) , (12) 



where. 



+, (t^u'^' + ^^T^A''^ (a^t/f - d,U; + g^f^U^U^) , (13) 

using T"-T'' = 2if'^'^T'^ in the case of SU(n) that is our interest in further discussion. 
Greek (Latin) indices run over 0,1,2,3 (1,2,3). Integrating out both sides in Eq. 
(12) over t (xi) for /i — i {/i — 0), and writing = [J'q, J), we then have, 

^U" - VUS = -gF j dt ( + F"^) = -^F j dx {J^, + F^) , (14) 
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which is satisfying the condition for irrotational fluid, i.e. the "vorticity" uj, 

a;" = V X U" = . (15) 

This condition is obtained as a non-trivial solution (U" 7^ 0) from the terms 
{djUf - diUf) = [V (V ■ U") - V^U"^]^. = [V x (V x U")]^. in Eq. (12), and 

djU^ — diUj = tjik (V X U")'^ in Eq. (13) for /x = j 7^ i, while it is vanishing for 
ji — i. Combining Eqs. (11) and (14), 

^ (7"v» + V (7» = -^F j dx ( Jp"o + ^0 ) ■ (16) 

At non-relativistic limit, 7 ~ 1 + l/2|vp and ~ 1 for a constant 0. These then 
yield, 

9v" 1_, „,2 



dt 2 



+ :tV|v1 



-^F^dx (Jp«o + Fo) 



(17) 

non— rel. 



up to 0(|vp) accuracy. Just to mention, one may also take a non-constant 0, for 
instance the hot fluid which takes 0non-rei.(^) ~ 1 + (5/2)(m/T) using large T/m 
expansion [9], inducing some additional terms in both sides of Eq. (17). Further, 
utilizing the vector identity |V |v|^ = (v • V)v -|- v x (V x v), finally we reach at. 



+ (v« • V) v'^ = -^F j dx ( Jp^o + F^) 



(18) 

non— rel. 



Obviously, Eq. (18) reproduces the classical EOM for irrotational fluid. We argue 
that Eq. (16) should be a general relativistic fluid equation. On the other hand, the 
"current force" JT"" is induced by the existing matters surrounded by and interacting 
with the fluid, while is induced by the fluid self-interaction and the interacting 
gauge flelds (A|^). Therefore, the lagrangian in Eq. (3) with the fluid fleld f/^ having 
a form of Eq. (11) should describe a general relativistic fluid system interacting with 
another gauge flelds and matters inside. Actually, the relative signs in the flrst and 
second terms in the left hand side of Eqs. (16), (17) and (18) rehes on the sign of U 
component in Eq. (11). This is the reason we should take the notation in Eq. (11) 
rather than the conventional one for relativistic velocity. 

To be more spcciflc, let us consider a magnetofluid uniflcation involving a non- 
Abelian SU(3) fluid contains particles of quarks and anti-quarks interacting with an 
electromagnetic field, i.e. the lagrangian with SU(3)f ® U(1)g gauge symmetry, 

C = iq^Q^-rnc^qq-^Sl^S^^^ -\f,,F^'' + g^J^^U^^ ^qJo^A^ . (19) 

from Eqs. (1), (4) and (5). Here qq is replaced with q denoting the quark charge, 
Q represents the quark (color) triplet, J^^ = QTp7^Q and Jq = Q7^Q, while Tp's 
belongs to the SU(3) Gell-Mann matrices. The model describes "macroscopically" 
non-Abelian fluid formed by dense gluon cloud surrounding the matters (quarks 
and anti-quarks) in an electromagnetic fleld. This flts the experimental clue from 
PHENIX Collaboration at the BNL-RHIC [2], far from being the weakly interacting 
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collisionless plasma the deconfined quark gluon matter behaves more hke a quark 
gluon fluid. This description is, however completely different with the hybrid mag- 
netofluid model [9, 10] where the fluid flows denote the quarks and anti-quarks 
interacting with the gluon fields. 

From Eq. (7) the zero-th component of fcrmion current, i.e. the current density, 
is J^pq — Jpo = Q}TpQ = Pp. Using Eq. (16) its relativistic dynamics obeys, 

^ (7 V» + V (7» = -QF j dx (p^ + F^) . (20) 

Meanwhile g-p should be associated with the fine structure constant of strong interac- 
tion = Anas with depends on the energy scale of physics under consideration 
[13]. According to the experimental results [3], the hot QGP is dense but seems 
to flow with tiny viscosity approximating the ideal fluid. This draws V • ~ 
affecting the term ff^d' {U^Uf) in Eq. (13), and a;" ~ satisfying Eq. (15). 
Moreover, as seen in the above lagrangian the gluon fluid itself does not feel the 
electromagnetic force, but the quarks and anti-quarks do. In contrary, the electro- 
magnetic force still contributes to the EOM in Eq. (20), but it is suppressed by a 
factor of q/ Qf oc ^Ja/as ~ O(10~^) since the quark charge is proportional to the 
electric charge e = ajA.'K. Therefore the electromagnetic force contributes about 
few percents, and might be neglected for a good accuracy. 

For completeness, we can derive the total energy due to fluid and electromagnetic 
fields in the above case. The lagrangian density for the fiuid and electromagnetic 
fields is, 

L = -\s;,S^'''' - \f^.F^^ + gpJ^^U^^ + qJG.A^^ (21) 

We can define Bp. = —Sq^ and ejikBp^ = Sj^ with i j denoting the electric- 
and magnetic-like fields of fiuid analogous to the usual electric and magnetic fields 
Ei = -Foi and ejikB'' = Fji. Then, S^^S"""" = -2 (|E^|^ - |B^|^) in analogy to 
Ff.uF^'' = -2 (|E|^ - |B|^). Due to the condition in Eq. (15), = |^f/"'"'U*xU^ 
while Ep = —qf § dx (pp + Fq ) from the EOM in Eq. (9). Hence, omitting the small 
electromagnetic current term we can easily deduce the Hamiltonian density, 

+ \ Ir^^U^ X Uf j + ^F7> (Pf + V" • J^) 

(22) 

In the simplest case of pure gluon fiuid, only the first and second terms will remain, 

+ ^ Ir^'U^ X U-^fj , (23) 

since ^ = in the absence of matter. Taking its non-relativistic limit, the total 
energy is simply H — J d^x7inon-rei.- For example having particular temperature 
dependent c/) would actually enable us to calculate any physical observables (free 
energy, entropy, etc) through the partition function density Z = Tr [exp(7i/T)]. 



n 




y"dx (p^p + FS) 



+ 2(|E|' + |Bn 



n 



dxF" 
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We have studied the magnetofluid unification using lagrangian approach and 
imposing the gauge principle in Abehan and non-Abehan cases. The model pro- 
vides new insights into macroscopic dynamics of relativistic hot fluids relevant for 
QGP. An example of SU(3)f<8) U(1)g model, in particular, may describe some of 
the puzzles posed by the data from the RHIC indicating that there is a very rapid 
thermalization in the collisions, after which a fluid with very low viscosity and large 
transport coefficients seems to be produced, and surprisingly the produced medium 
seems to be strongly interacting [2]. Recently, the ALICE experiment at the LHC is 
also challenging to resolve this puzzle [14] . It is also argued that the formalism leads 
to different physical consequences than the other ones derived from the hybrid mag- 
netofluid model. However, all consequences here arise naturally from flrst principle 
without any flne tuning. The results should lead to some interesting phenomenon 
either in QGP or fluid based cosmology. Such studies, including simulations based 
on numerical computation for particular cases, are in progress. 

We greatly appreciate fruitful discussion with N. Riveli, H.B. Hartanto and A. 
Oxalion throughout the work. This work is partially funded by the Indonesia Min- 
istry of Research and Technology and the Riset Kompetitif LIPI in flscal year 2007 
under Contract no. 11.04/SK/KPPI/II/2007. 
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